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ABSTRACT 

Rates of coke formation during naphtha pyrolysis have 
been investigated in a jet-stirred reactor both for thiophene- 
free naphtha and thiophene containing naphtha on Inconel 600 
surfaces at atmospheric pressure in the temperature range of 
1073-1103 K, The effect of temperature# space time# weight 
ratio of steam to naphtha# aromatic and thiophene content of 
feed and material of construction on the rates of coke deposition 
was measured by periodically weighing a snail metalplate 
suspended into the centre of the reactor. 

Coking rates increased as the temperature# space time 
and aromatic content of the feed was increased whereas it 
decreased with an increase in the weight ratio of steam to 
naphtha and the thiophene content of the feed. Addition of 
thiophene significantly reduced the aromatic yields. Rate of 
coke formation could be modelled using a power law form. 



CHAPTER 1 


INTRODUCTION 

Pyrolysis of hydrocarbons, such as naphtha, ethane and 
propane is of great importance for the production of olefins 
such as ethylene, propylene and butadiene which are the build~ 
ing blocks of the petrochemical industry. The preferred feed- 
stock varies from country to country. In the U.S,, it is 
mainly a mixture of ethane and propane, whereas in Japan, 

Europe and India, the preferred feedstock is naphtha. Naphtha 
is a complex mixture of hydrocarbons and is normally characterised 
with respect to its specific gravity, carbon to hydrogen ratio, 
ASTM boiling range and PONA analysis, which provides the 
distribution of paraffins, olefins, naphthenes and aromatics. 

A major problem associated with any hydrocarbon pyrolysis is 
the undesirable formation of coke which deposits on the walls 
of the pyrolysis coil. Coke is a complex mixture of carbon 
and hydrogen and can be formed either from the reactants and/or 
products. The deposition of coke with time gradually 

increases the pressure drop across the tubular reactor and 
also reduces the overall heat transfer coefficient across the 
tube wall due to the extra thermal resistance of the deposited 
coke layer. Hence, the thicker the coke layer, the higher 
the tube wall temperature which must be applied to achieve 
the sane fluid temperature and degree of feedstock conversion. 
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The reactor has to be shut down when the tube wall temperature 
reaches the limiting value. Depending upon the feedstock and 
operating conditions, the run lengths are generally limited 
from 40 to 60 days. The reactor is then decoked by burning 
with a steam-air mixture. Frequent decoking operation results 
in loss of production, affect the coil life and utility costs 
are increased. Coke is also corrosive and can cause metallur- 
gical failure of a reactor. 

To design and optimise naphtha crackers, reliable models 
for the main naphtha pyrolysis as well as models for the 
Coking kinetics during pyrolysis are necessary. Not much 
information is available on modelling of naphtha pyrolysis or 
coking kinetics. Quantitative models are not available, but 
only the qualitative effect of various operating variables is 
known. 

During the last few years, investigations have been made 
on coke formation during the pyrolysis of various hydrocarbons 
such as ethane (Albright and McConnell, 1978; Sundaram et al,, 
1981), propane (Sundaram and Froment, 1979; Trimm and Turner, 
1981a), n— hexane (Pramanik and Kunzru, 1985), n-octane (Shah 
et al,, 1976) and naphtha (Newsome and I«ftin, 1979; Kumar and 
Kunzru, 1985), These studies have shown that coke deposition 
depends on several factors, such as, the aromatic and sulphur 
content of the feedstock, hydrocarbon partial pressure, tempera- 
ture, conversion and the material of construction of the 


reactor 
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In an earlier study, (Kumar and Kunzru, 1985), coke 
deposition in naphtha pyrolysis was studied in a jet stirred 
reactor and the rate of coke foirmation was modelled by an 
approximately second order reaction involving the aromatics. 
However, the model was not checked either for different 
aromatic concentrations of the feedstock or for different 
materials of construction of the reactor. Moreover, the 
effect of sulfur content of the feed on the rate of coke 
formation was not investigated. It has been observed that the 
continuous dosage of sulfur containing compounds such as 
thiophene, together with heptane feed, influences both the 
rate of pyrolysis as well as the rate of coking (Bajus et al,, 
1981), However, it is not clear whether the sulfur caused 
passivation of the reactor walls or it affected the coke 
precursors in the gas phase. 

The objective of this study was to investigate the effect 
of the material of construction of the reactor and the aromatic 
and sulfur content of the feed stock on the rates of coke 
formation during naphtha pyrolysis. In addition, the effect 
of temperature, inlet hydrocarbon partial pressure and space 
time on the rate of coke formation was investigated. A further 
objective was to model the coking kinetics during naphtha 
pyrolysis. 



CHAPTER 2 


LITERATURE REVIEW 


Product distribution in the pyrolysis of various hydro- 
carbons including naphtha is well established and considerable 
data is available, but very meagre data is available on the 
rates of coke formation during the pyrolysis of naphtha or 
even other pure hydrocarbons. 

Various investigators have studied the pyrolysis of 
methane (Eisenberg and Bliss, 1967), ethane (Silcocks, 1956/ 
Snow et al., 1959), propane (Volkan and April, 1977), n-butane 
(Steacia and Puddington, 1938; Sandler and Chung, 1961), 
n-hexane (Frey and Hepp, 1933; Ebert et al,, 1983), n-heptane 
(Bajus and Vesely, 1979), n-nonane (lOanzru et al,, 1972), 
2-pentene (Kunzru et al,, 1973), naphtha (Van Damme et al,, 
1981; Kumar and Kunzru, 1985) and have reported the kinetics 
and product distribution. 

Data on coke formation during pyrolysis is very limited, 
Hirt and Palmer (1967) studied the coke formation from methane 
pyrolysis between 1163 and 1373 K and obtained an activation 
energy of 432 MJ/kmol, Johnson and Anderson (1972) studied 
the coke formation from acetylene at temperature ranging from 
773 to 1273 K and observed both coke and polyma: as the 
products. They suggested that coke and polymer were formed 
by two competing parallel reactions in the gas phase, Albright 
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and McConnell (1978) Investigated the rate of coke formation 
in ethane pyrolysis and found that the coking phenomena was 
significantly affected by the material of construction of the 
pyrolysis tube and also on the pretreatment given to the 
reactor. An activation energy of coke formation of 250 to 
290 MJ/kmol was reported, 

Kinney and Del Bel (1954) studied the pyrolysis of unsatu— 
rated aromatic hydrocarbons and suggested that coke is formed 
from a stepwise condensation of the aromatic rings. Similar ly# 
Virk et al, (1974) have proposed that coke formation occurs 
from aromatic compounds, mainly from benzene, Shah et al, (1976) 
studied the coke fortnation in the thermal cracking of n-octane 
between 1023 to 1073 K. The coking rate was high during the 
Initial stages of the experiment and then reduced to an 
asymptotic value. They also observed that the amount of coke 
formed increased with an increase in the surface to volume 
ratio of the reactor. Further air (or oxygen) treatment activate 
the reactor surface to induce more coking whereas H 2 S treatment 
passified it. The reduction in coke formation during hydro- 
carbon pyrolysis by adding sulfur compounds to the feed is 
well documented (Bajus et al., 1981), Newscme and leftin (1979) 
investigated the rates of coke formation in naphtha pyrolysis 
in an annular tubular reactor. The coking rates were found 

to depend on the run time, inlet steam and hydrocarbon partial 
» 

pressure and temperature. The spjecific rate of coking was 
independent of the surf ace-to-volume ratio of the reactor. 
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All the above mentioned studies were conducted in 
tubular flow reactors and the total coke formed measured 
by either weighing the reactor before and after the experiment 
or by burning the coke and measuring the effluent gaseous 
amounts. Since coke is not deposited uniformly along the 
reactor (Froment and Bischoff, 1961), this result is only 
average coking rate and precludes a kinetic analysis. To 
overcome this problem, Sundaram and Eroment (1979) developed 
a jet-stirred completely mixed reactor which allowed the 
kinetics of the main pyrolysis reaction and kinetics of 
Coking to be determined simultaneously. They measured the 
amount of coke deposition continuously, by suspending a small 
cylinder of the same material as the reactor, into the reactor 
from an electrobalance. For propane pyrolysis, the coking 
rates were best represented by a model which assumed that 
coke was formed from propylene by a first order reaction, 
whereas the coke formation in ethane pyrolysis was modelled 
as a first order reaction from the C^+ hydrocarbons 
(Sundaram et al, ,1981 ) ,Trimm and Turner (198la) also studied 
the coke formation in propane pyrolysis in a jet-stirred 
reactor and found that once the reactor surface was covered 
with carbon, the rate of coke deposition was independent of 
the nature of the surface at all temperatures, Pramanik and 
Kunzru (1985) investigated the pyrolysis of n-hexane and 
kinetics of coking in a jet -stir red reactor In the temperature 
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range of 993—1083 K, The experimental coking rates could 
be adequately fitted by a model in which ethylene was the 
coke forming species. Kumar and Kunzru (1985) also used a 
jet-stirred reactor to study the coke deposition during 
naphtha pyrolysis and the rate of coke formation was modelled V- 
by an approximately second order reaction involving the 
aromatics. 

The reactor surface significantly affects the coke 
formation during pyrolysis of hydrocarbons. Various investi- 
gators have shown that metals present in the reactor surface 
activate the coke formation. The presence of metals in coke 
during ethylene cracking has been reported by Albright and 
coworkers (1978), who have demonstrated that such metals can 
catalyze the formation of additional coke. Research of Albright 
and co-workers (1979) on the pyrolysis of hydrocarbons over 
Incoloy 800 surfaces shows iron to be the predominant metal in 
the coke. 

Since it is well established that the walls of the reactor 
catalyze the coke formation, most of the methods used to inhibit 
the coking use some means to passivate the reactor walls. 
Presulfiding the reactor has been reported to reduce coke 
formation during the pyrolysis of ethane (Alright and McConnel, 
1978), propane (Trimm and Turner, 1981b) and octane (Shah et al,, 
1976). Bajus et al, (1981; 1983) studied the effect of thio- 
phene, dibenzyl sulfide and dibenzyl disulfide on reaction 
kinetics and coking of steam cracking of heptane in stainless 
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steel tubular reactors with relatively large surface areas. 

The deposition of coke was inhibited considerably due to 
the formation of metal sulfide which passivated the reactor 
surface, Bajus et al, (1980) studied the effect of sulfur 
compounds on straight run naphtha and proposed that during 
the cracking, intermediate HS* radicals could be formed from 
sulfur compounds, such as mercaptans. These radicals might 
be partly converted to hydrogen sulfide by reaction with 
hydrocarbons and could also react with the reactor walls 
to form metal sulfides. This sulfide film inhibits the 
secondary reaction which leads to the formation of coke, 

Depeyre et al. (1985) studied the influence of sulfur compounds 
on n-nonane in a quartz reactor and found that the addition 
of sulfur compounds to n— nonane steam cracking catalyses the 
rate of n-nonane decomposition and inhibits secondary reactions 
leading to the formation of C 2 H^, aromatic compounds. 

Yields of these secondary products decreased as the weight 
percentage of sulfur compounds in the feed was increased 
and sulfur also altered the rate of formation of coke. 

From the above literature survey, it is clear that not 
much information is available for the rate of coke formation 
in naphtha pyrolysis and very few quantitative models have 
been published. Moreover, the effect of sulfur compounds on 
the rate of coke formation during steam cracking of naphtha 
is not well established. 



CHAPTER 3 


EXPERIMENTAL APPARATUS AND PROCEDURE 


3 *1 Exper Im ent a 1 Appar at u s t 

The experiments were conducted in a jet— stirred reactor, 

A schematic diagram of the equipment is shown in Fig, 1, 

Liquid naphtha and water from separate graduated burettes 
were pumped by a double acting reciprocating pump which could 
be adjusted for flow rates in the range of 0,2 to 30 cc/min. 
Steam, which was used as an inert was generated in a vaporiser 
(1,5 Kw) and mixed with naphtha. The mixture of hydrocarbon 
and steam was passed through another preheater (1,5 Kw ) where 
the mixture was heated to approximately 773 K, The heated 
mixture was then led into the preheater which was just above 
the reactor. The product gases from the reactor exited 
through a central outlet tube which extended above the preheater 
and was constructed from 25 mm 0,D, SS 316 tubing. 

The preheater reactor assembly was such that good mixing 
could be achieved in the reactor. It was made from a 25 mm. 

I.D. SS 316 tubing and the reactor length and volume was 2,3 an 
and 11.2 cm^ respectively. The wall separating the reactor and 
preheater contained thirty two one mm holes, distributed 
uniformly, which acted as nozzles and ensured good mixing 
inside the reactor. To provide uniform distribution of the 
Inlet mixture into the reactor , the preheater was packed with 






Fig.t Schematic diagram of the experimental setup 
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crushed ceramic beads. The preheater-reactor assembly was 
heated in a controlled electric furnace. The temperature 
of the reactor was measured and controlled by a temperature 
indicator-controller (Model CT 806, Century Instrument) 
which was connected to a chromel— alumel thermocouple placed 
in a thermowell which extended into the centre of the reactor. 
The reactor temperature could be controlled to within + 2®C 
of the set value. There was a large temperature gradient 
in the preheater zone, but, in all the runs, the axial 
temperature gradient in the main reactor was negligible. 

In this set up, both the deposition of coke and product 
distribution could be measured. For the measurement of coke, 
a small rectangular metal plate was suspended by means of a 
wire into the central portion of the reactor through the 
central outlet tube. The plate was weighed on an electro- 
balance and resuspended at the centre of the reactor. This 
was done at periodic intervals during the course of a run. 

To prevent reaction products from entering the outlet tube, 
a small flow of nitrogen was passed from the top through the 
side exit of the outlet tube. The reaction mixture was quenched 
by passing through two ice-water cooled condensers placed in 
series. The liquid product was collected and the non-conden- 
sables passed through a sampling valve and vented, 

3,2 Ebcperimental Procedure ; 

A run was initiated by first heating the reactor to a 
temperature 1 0-1 5 *G higher than the desired temperature. 
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After the temperature had stabilised, nitrogen gas was 
passed from the top. The reactor was then flushed with 
steam for 45 minutes to ensure an inert atmosphere in the 
set-up. The water flow rate was then fixed to the desired 
value and naphtha feed started. Since pyrolysis is endo- 
thermic, the reactor temperature dropped by 10— IS^C once 
the reaction commenced. Depending upon whether the product 
distribution or coke deposition was being measured, the 
procedure was modified. 

For measuring the coke deposition, after the tanperature 
had restabilised at the new value, a small plate was suspended 
by, means of a wire into the central portion of the reactor 
through the central outlet tube. The suspended plate was 
removed from the reactor at thirty minute intervals, weighed 
on an electrobalance and resuspended. The accuracy of the 
balance was 0.01 mg. To prevent any oxidation of the deposited 
coke, the plate was pulled up slowly through the outlet 
tube so that it cooled to 473-523 K in the relatively cooler 
portion of the outlet tube, A run was continued till the rate 
of coke formation became constant. The total duration of the 
run varied from 2-3 hours. Attempts to measure the coke 
deposition by suspending the plate directly from an electro- 
balance were unsuccessful because of vibrational problems 
during the passage of reactants and products, .IXiring a run, 
the total flow rate of the non-condensables was measured. 
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At the ccxnpletion of the run, liquid products were collected 
and the volume and density measured. The aromatic content 
of the liquid products was determined by gas chromatographic 
analysis. Subsequently, the reactor was flushed with steam 
and the deposited coke on the plate and reactor was burned 
by passing heated air through the reactor for approximately 
30 minutes. The plates were used for the next run after 
removing the metal oxides by immersing the plate in a solution 
of acids (3% H2S0^ + 3% HCl) and then cleaning with emery 
paper . 

For measuring the product distribution, no plate was 
suspended and the top outlet tube was closed with a cork. 

The total gas and liquid flow rate was measured and the gas 
from the sampling valve analyzed on a gas chromatograph. 

The analysis of the liquid products was done as mentioned 
above, 

3,3 Liquid Analysis : 

Due to the complexity of the liquid products and the 
nonavailability of a suitable capillary column, detailed 
analysis of the liquid fraction was not possible. The 
aromatic content of the liquid product was determined on a 
3 mm I,D., 3m long DEGS (diethylene glycol succinate) column. 
The calibration data for liquid analysis is ^own in 
Appendix A, Details of the chromatographic conditions are 
given in Table 1, 
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3 #4 Gas Analysis I 

The gaseous products formed in the pyrolysis were 
analysed on gas chromatograph using a thermal conductivity 
detector, A 3mm I,D. , 3 m long Duropak column operated at 
42 ®C was used to analyse the and gases. Hydrogen, 
methane, ethane and ethylene could not be completely 
separated on this column and were analysed on a 3 mm I.D, , 

2 ra long poropak Q column also operated at 42 ®G, The cali— 
tration data for the gas analysis on the different columns 
is shewn in Appendix A, Etetails of the chromatographic 
conditions are given in Table 1, 



Table 1; Details of the Chromatographic Condition 


15 




% 




\D 


m 


X 


C 0) 


03 

U 






C (U 


% 

% % 


o 

rH 

K K 

X 

0) 0) 

C C <0 

6 0^ 

0 0) 

u 

(U 0) c 

N D Cl) 

U W 

CO ^ 

U U 

^ X 

03 04 

C rHI rH 

0) O >1 


D 

H 

PQ £• X 

4J 




C jj 




0) c 




p (U 




B M c 

M E 

o 

o 

O 

-H 0 

a» 

(y\ 

(T\ 

vO 




rate 

i 

X 10 

oo 

\D 

0- 

04 

CO 

vO 

vD 

IS w 

CD 

vO 

vD 

O \ 

CN 

If) 

vO 

r-l iw ro 

• 

• 

• 

1^0 H 

o 

o 

o 

o 

4J ^ 

O 

o 

o 

n • 

• 

• 

• 

0) 

CO 

CO 

CO 

E 



tH 

C 0) 

H 4J 

to 

CO 


0 

4^ ^ 

c 

o 

O 

V • 

• 

e 

• 

35 p, - 

ro 

CO 

CO 

4-> E 



T— < 

0 55 

CO 

CO 


W 4J 








• 




C 04 

O 

o 

G 

0) E « 

• 

• 

• 

> 0) 

m 

in 

CO 

O -M 

t-i 

rH 

vO 


ro 

to 

CO 


X 

X 

X 

n 




£ (U 

• 

• 

• 

§ N 

Q 

P 

P 

rH ^ 

• 

' • 

« 

0 to 

H 

H 

H e 

o 

E 

£ 

E ro 


E 

E 


E ro 

E ro 

E 


CO 

CO 

CO 




A 




C 

cn 


a 

O A 

CC ' 


m 

14 

' E -H , 


% 

CD O 

3 X 


X 

O to 

iH a 

IB 


go 

o to 

A 

A 

P E 

o A 

0 

O 


i4 

u 

^ U 


3 

Q 

in jc 


P 

A 

04 o 








CHAPTER 4 


RESULTS AND DISCUSSION 


To study the coke formation during naphtha pyrolysiS/ 
experiments were conducted at atmospheric paressure in the 
temperature range 1073 to 1103 K, The naphtha was obtained 
from Bharat Refineries^ Bombay ^ which processes Arabian crude 
and the feed characteristics are given in Table 2, Steam was 
used as an inert diluent and the weight ratio of steam to 
naphtha^ 5 t was varied from 0,47 to 0,7 and the space time 
from 0,2 to 0,49 seconds. The effect of aromatic and sulfur 
concentration of naphtha feed on the rate of coking was studied 
by adding benzene or thiophene to the naphtha feed. The 
benzene concentration was varied from 8,1 to 21,07 mole % 
and the thiophene content from O to 1,2 wt,%. The rate of 

coking was investigated on three metal surfaces viz, SS 304, 

SS 316 and Inconel 600, For all the runs, a material balance 
was made and the water balance was 98-99% and the hydrocarbon 
balance was 96 — 98%, 

4,1 Effect of Material of Construction t 

The effect of material of construction was investigated 
on SS 304, SS 316 and Inconel 600, For all the runs, the 
temperature, 5 end space time, ^ / were maintained at 1083 K, 

0,47 and 0,49 s, respectively. The variation in the rate of 
coke formation with run time is shown in Fig, 2, For all the 
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Table 2: Characteristics of the Naphtha Feed 


Properties ; 

(i) Density (298 K) 0.701 gm/tnl 

(ii) ASTM Distillation *C 

IBP 47 

5% 59 

10% 63 

50% 87 

90% 119 

95% 129 

FBP 156 

(iii) PONA wt.% 

Normal paraffins 33,7 

Total paraffins 73,6 

Naphthenes 17,4 

Benzene + Toluene + Xylene 8,5 

(iv) Total sulfur' 0,05 


Average molecular weight 


91 
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runs, rates of coke deposition are initially high and then 
decrease to a constant value with progress of the reaction, 

A plausible explanation for this variation is that initially 
the bare metal surface catalyzes the coke formation and as 
the coke gradually deposits on the surface, the metal 
catalytic activity is progressively reduced. The asymptotic 
coking rate period, then, corresponds to coke deposition on 
the metal surface completely covered with coke. Similar 
trends for rates of coke formation have been reported during 
the pyrolysis of n-octane (Shah et al,, 1976), naphtha (Kumar 
and Kunzru, 1985), propane (Trimm and Turner, 19813) and n-hexane 
(Pramanik and Kunzru, 1985), As shown in this figure, the 
coking rates were highest on SS 304 and lowest on Inconel 600, 
This indicates that coking rate is dependent on the material 
of construction of the reactor, Mandal (1984) also reported 
coking rates during n-hocane pyrolysis to be lower on Inconel 
600 than on SS 304, The dependency of coking rates on the 
nature of the metal surface has also been reported by Brown 
and Albright (1976), Albright and McConnel (1979), Contrary 
to these findings Trimm andTurner ( 1981a) reported that although 
the initial coking rates were different, the asymptotic coking 
rates were identical for different surfaces. Since, in this 
study, the coking rates were significantly lower on Inconel 600, 
subsequent investigations were done only on Inconel 600 


surface 
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4,1,1 Product Distribution 

For some runs on Inconel 600, the major gaseous products 
were also determined. As expected, the major products were 
ethylene, propylene, methane, together with small amounts of 
hydrogen, butane and ethane. The aromatics constituted a 
major portion of the organic liquid products. The product 
distribution at 1083 K and a space time of 0,49 s is shewn in 
Table 3, These yields compare favourably with those reported 
in the literature (Kumar and Kunzru, 1985). To check the 
reproducibility of the product distribution and coke deposition 
measurements, an identical run was repeated three times. The 
rates of coking were reproducible to within + 1% and the 
product distribution to within + 2%, 

4,2 Effect of Decoking : 

To study the effect of metal oxides formed during 
decoking on the rate of coke formation, an identical run was 
made first on an Inconel surface from which the oxides had 
been removed by treating with 3 % H2S0^ + 3 % HCi solution and 
then with the oxide layer formed during decoking still present. 
The variation of the rate of coke formation with run time 
for the two cases is shown in Fig, 3, As shown in this 
figure, the presence of metal oxides increases the rate of 
Coke formation. The effect of the metal oxide was reproducible 
as checked by again measuring the rate of coke fcarmation on 
the decoked (but not cleaned) surface. Since the rates were 
identical, these points have not been ^ewn in Fig, 3, 
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Table 3: Product Yields During Naphtha Pyrolysis 


Temperature = 1083 K; 

6 0,47 Kg Steam/Kg naphtha 

T = 0.49 s 



hydrogen 1,20 
methane 15,92 
ethylene 25,32 
ethane 2,36 
propylene 7,32 
n-butane 0,05 
other gases 28,83 

Benzene + Toluene + Xylene 14,7 
high^ boiling liquid 2,85 
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There Is scxne controversy regarding the effect of the 
metal oxides on the rate of coke formation. Brown et al,, 

(1982) mention that commercial ethylene plants normally rely 
on preoxidation during the steam/air decoking operation to 
produce a protective surface oxide film. On the other hand 
Brcwn and Albright (1976) reported that metal oxides promote 
secondary and undesired reactions leading to coke, Bajus et 
al, (1981) found that in order to achieve reproducible results, 
the metal oxides formed during decoking had to be removed by 
the action of solution of acids prior to each run, Albright 
and Yu (1978) reported that metal oxides on the surface are 
likely to favour the production of active coke. 

In this study, the results were reproducible both for 
cleaned Inconel 600 and Inconel 600 with metal oxide present. 
Since the rates of coke formation were lower when the oxides 
had been removed by acids, all the data were taken on surfaces 
from which oxides had been removed, 

4,3 Effect of Aromatic Concentration : 

It is well known that the nature of the feedstock affects 
the coke deposition (Ranzi et al, 1985), Although the role of 
benzene and other aromatics in increasing the rate of coke 
formation during pyrolysis is well documented (Virk et al,, 

1974? Kumar and Kunzru, 1985), not much quantitative information 
is available regarding the effect of aromatic concentration of 
naphtha on coking rates. To study the effect of aromatics on 
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the rate of coke formation* The aromatic concentration of the 
naphtha feed was varied from 8,1 to 21,07 mole% by adding 
benzene to the naphtha. The effect of inlet aromatic 
concentration on the rate of coke formation is shown in Fig, 4, 
Rates of coke formation increased with an increase in the 
inlet benzene concentration. For instance, as the inlet aro- 
matic concentration was increased from 8,1 to 21,07 mole%, 
the asymptotic coking rate increased by 220%, Most probably, 
the increase in coking rate is due to the formation of unsatu- 
rated compounds with the elimination of hydrogen in the 
condensed ring (Kinney and Del, 1959; Virk et al,,‘ 1974), It 
should be mentioned that at these conditions benzene does not 
pyrolyse. 

To study the effect of the nature of the aromatic, one 
run was made with 11,49 mole% toluene, instead of benzene, 
added to the feed. For the same concentration the effect of 
toluene and benzene on the rate of coke formation was identical, ; 

4.4 Effect of Process Variables ; 

The rates of coke formation varied with reactor temperature^ 
inlet naphtha partial pressure and space time. To study the 
effect of temperature on coking rates the reactor temperature 
was varied from 1073 K to 1103 K which is the usual range for 
naphtha pyrolysis. Since the variation of coking rates with 
run time was similar to that shown in Fig, 2, only asymptotic 




Rate of coking X 10^, kg/(m^)(h) 
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coking rates are being reported. The detailed variation of 
coking rates with time for the different runs is given in 
Appendix B, The variation of the asymptotic coking rate, r , 
with temperature and space time is shown in Fig. 5. Coking 
rates increase with increasing temperature and space time. 

This trend is similar to that reported for the pyrolysis of 
propane (Trimm and Turner, I98le)^ n-hexane (Pramanik and Kunzru, ! 
1985) and naphtha (Kumar and Kunzru, 1985), Since coke is 
formed from secondary reactions and once formed does not react 
further, rate of coking increases with temperature and space . 
time. 

The effect of inlet steam to naphtha ratio on the rate of 
coking at a constant temperature of 1083 K is shown in Fig, 6. 
With decreasing partial pressure of naphtha, the concentration 
of coke forming species is decreased, thus reducing the rate 
of coke formation. At these temperatures, the coke-steam 
reaction is not expected to be significant (Biba et al, 1978) and 
would not contribute to the reduction in rate of coking with 
increasing partial pressure of steam. These trends are similar 
to that reported by Kumar and Kunzru (1985) who used a 5 of 
0.71 to 1.43. 

4.5 Effect of Thiophene Content of Feed : 

The effect of sulfur content of the feed was studied by 
adding thiojiierte to the najiitha. The addition of various 
quantities of thiophene was so arranged that the wt .% of sulfur 
expressed as total weight of feed, varied from 0.05 to 0.5%, 









29 


It may be noted that the naphtha sample contains 0,05 weight % 
sulphur. 

The addition of thiophene affected both the product yields 
as well as the rate of coke formation. The rate of coke 
formation decreased with an increase in the sulfur content 
of the feed and the variation with run time for different 
thiophene concentrations is shown in Fig, 7, whereas Fig, 8 
shows the variation of the asymptotic coking rate with the 
thiophene content of the feed. As shown in Fig. 8, till 0,4 
wt,% thiophene, the asynptotic coking rate was inversely 
proportional to the thiophene content of naphtha whereas at 
higher thiophene concentrations the inhibition in coking rate 
was less. 

To check whether the addition of thiophene had an effect 
on the primary pyrolysis reaction, the total gas and liquid 
yields for different thiophene concentrations were measured. 

In addition, the aromatic content of the liquid product and 
the major gaseous product were also analysed. The addition of 
thiophene to naphtha was found to affect both the total gas 
and liquid yields. As shewn in Fig, 9, the total gaseous 
yield increased whereas the total organic yields decreased 
with an increase in the thiophene content of the feed. The 
effect on the gaseous yields was highest at lower thiophene 
concentrations whereas between 0,4 — 1,2 wt.% thiophene, the 
gas yields were directly prop>ortional to the thiophene concen- 
tration, For the rarge of thiophene concentration investigated. 
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Table 4: Product Distribution During Pyrolysis of Naphtha Containing Thiophene 
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the liquid yields were inversely proportional to the thiophene 
concentration of the feed. 

Although the composition of the gaseous products did not 
change due to the addition of thiophene in the feed, the yields 
of the gaseous components, expressed as wt,% of naphtha feed, 
was higher because the total gaseous yields were higher. The 
effect of thiophene addition on the product yields is shown 
in Table 4, Addition of thiophene significantly reduced the 
aromatic content of the liquid product. Addition of 0.13 % 
thiophene in feed reduce the aromatic yield from 14,7 to 
5,78 wt,%. On further addition of thiophene the decrease in 
the aromatic yield was less significant. As the thiophene 
concentration was increased from 0,13 to 1,2 wt,% the aromatic 
yield decreased from 5,78 to only 4,18 wt,% , It should be 
noted that the aromatic yield with thiophene is less than the 
aromatic content of the thiophene free naphtha feed, which 
indicates that at this temperature thiophene is reacting with 
the aromatics. Aromatics are very refractory and in general, 
do not crack during naphtha pyrolysis (VanDamme et al,, 1981). 
However with thiophene present in the reaction mixture, the 
aromatics can react with S. and .SH free radical produced due 
to thiophene decomposition (Glass and Reid, 1929), 

The sulfur present in the feed can reduce the rate of 
coking in two ways. It can either react with th6 metal surface 
to form metal sulfides thus passivating the reactor walls or it 
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can take part In the homogeneous gas phase free radical 
reactions {Bajus et al., 1981). To check whether sulfur 
was taking part in the heterogeneous and/or homogeneous 
reactions one run was conducted with 1,.2 wt,% thiophene in 
the feed and after the asymptotic coking rate had been attained, 
the reactor was flushed with steam and thiophene-free naphtha 
was passed without any intermediate decoking. It was found 
that the asymptotic coking rate was identical to that 
measured when a fresh Inconel 600 surface was used. This, 
together with the change in product yields, shows that thiophene 
is definitely taking part in the homogeneous reactions. If 
the reduction in rate is solely due to the presence of metal 
sulfides then the rate for the two cases should have been the 
same, because the passivating metal sulfide layer had not 
been removed by oxidation during decoking. The effect of 
sulfur In passivating the metal wall is well known (Bajus 
et al,, 1980; Trimm and Turner, 1981b) and as such no attempt 
was made to confirm the presence of metal sulfides. 

These results are in good agreement with the limited 
information available regarding the effect of sulfur ccxnpounds 
on the rate of coking, Bajus et al, (1981) studied the 
effect of adding thiophene to the raffinate feed and reported 
that addition of thiophene affected the coke formation as well 
as the pyrolysis. They also observed that till 0,5% thiophene, 
the rate of coking was Inversely proportional to the thiophene 
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content and production of organic liquid products decreased 
with an Increase In the thloj^ene content of the feed. During 
n-heptane pyrolysis, Bajus et al, (1981, 1983) found that the 
addition of thiophene, dlbenzyle sulfide or dlbenzy3^ disulfide 
decreased the rate of coking. At low sulfur concentration 
of the feed, the rate of n-heptane pyrolysis was Increased 
whereas a further Increase In the sulfur concentration 
decreased the rate of pyrolysis, Depeyre et al, (1985) reported? 
the addition of sulfur compounds to n-nonane steam cracking 
catalyzed the rate of n-nonane decomposition and Inhibited 
secondary reactions leading to the formation of <^ 4 ^^ and 

aromatic compounds. Yields of these secondary products decreased: 
as the weight percentage of sulfur compounds In the feed was 
Increased, Sulfur reduced the rate of formation of coke. 

Comparison of kinetics, selectivity of decomposition and 
rate of coking of naphtha leads to the finding that thiophene 
influences both the radical conversion in the homogeneous 
phase and the surface reactions in the reactor. Under the 
given experimental conditions, reactions can occur between 
thiophene molecules themselves and between thiophene and the 
highly reactive radicals resulting from the decomposition 
of hydrocarbons, e,g, H, and R, 

* 

► 
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Wynberg and Bantjes (1959) pyrolyzed thiophene at 1073 - 
1123 K, Among the reaction products, bithienyl as well as 
benzothiophene and phenylthiophene was found. The formation 
of these products was explained by Diels-Alder addition of 
thiophene itself according to the reactions (2) and (3). 



Fields and Meyerson (1967) ascribe a very important role of 
thiophene, which can be formed by intramolecular dehydrogenation 
of thioj^ene, as a reactive intermediate product. Reaction 
with a thiophene molecule could give benzothiophene according 
to reaction (4), 



From the proposed mechanisms for the thermal conversion of 
thiophene, the release of sulfur (reactions 1, 3, 4) and hydrogen 
sulfide (reaction 2) is characteristic. Sulfur eliminated from 
thiophene can enter into reactions in the gaseous state. It 
can also react with the metal surface of the reactor. From the 
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possible reactions / sulfur can initiate propagation according 
to reaction 5 with the formation of ,SH radical. 


C + S 

n "2n+2 




(5) 


Radical *SH from reactions 2 and 5 can also cause propagation 
according to reaction 6, termination according to reaction 7 
or react with the metal surface to form metal sulfides 
according to reaction 8. 


^n «2n+2 


+ , *SH 



«2n+l 


+ H2S 


( 6 ) 


*SH + 'h H 2 S (7) 

*SH + metal ■ »- metal sulfides + H* (8) > 


Reactions such as (5) and (6) would increase the rate of 
pyrolysis of naphtha resulting in higher gaseous and lower 
liquid yields. 


4,6 Modelling of Coke Formation : 

Coke formation in the pyrolysis of naphtha is a complex 
phenomena due to the various possible coke forming free radical 
reactions and the exact mechanism is not clear, 

4,6,1 Modelling of Coke Formation During Pyrolysis of 
Thiophene»Free-Naphtha : 

Kumar and Kunzru (1985) postulated simple models involving 
either the reactant and/or products to model coke formation 
during naphtha pyrolysis. The feed was thiophene-free but the 
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water used contained 50 ppm CS 2 . On SS 304 surface, the 
asymptotic coking rate could be represented by the following 
expression 

= 1.68 X lO^^ . exp ( ^.1.98 

where r^ is the asymptotic rate of coking in (kg coke )/(hr ) (m^ ) 

and is the concentration of the aromatics in the reactor in 
3 

kg mole/m , For this study also the applicability of relating 
checked. The model postulated was 

exp (. 1^) (2) 

Since the experiment was conducted in a mixed reactor, the 
concentration of aromatics in the exit was the same as that 
inside the reactor. For calculating the molecular weight 

of the gaseous product was taken to be 30 (Kumar and Kunzxu, 
1985) and the molecular weight of the aromatic free liquid 
product to be 91. Since the total volumetric flow rate of 
the gaseous product and the mass flow rate of the liquid product 
was experimentally measured, could be calculated as, 

^ _ moles of aromatics 

A ” moles of aromatics + molesof aromatic 
free product + moles of water + moles 
of gaseous product 

For all the runs, the molar flow rate of the aromatic - free 
product constituted only from 2.0 to 4,0% of the total 
molar flow rates. As shown in Fig. 10, it was possible to 
relate r^ for all the runs with C^. It should be emphasized 
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that the points shown in Fig, 10 were taken at various 
temperatures, space times, 6 , and inlet benzene concentration 
of the feed. 

The initial estimates of k , E and n were obtained 

o 

graphically using linearization. To estimate the order of the 

reaction In r^ vs In was plotted for various temperatures 

and the plot is shewn in Fig, 10, The reaction order for coking 

varies from 1,0 to 1,4 and the average order for the rate was 

1,29, However, for modelling purposes the value of n was 

taken to be a constant. The rate constants at different 

temperatures were obtained from the slopes of the plot between 
1,29 

r^ and C^' , as shewn in Fig, 11, The value of the activation 

energy using linear analysis was obtained by an Arrhenius plot. 
Fig, 12, of Ink vs 1/T, and the value was found to be 58,2 kcal/ 
gmole. 

The parameters of expression (2) i,e, E, and n were 
estimated by nonlinear optimization using data for 23 runs and 
the tabulated values of r^, T and are shewn in Table 5, 

The objective function was chosen to be 

Z^C^^a^cal*” ^^a^expt,] 
where k^ exp (- |f> 

The objective function was minimized by the appropriate choice 
of frequency factor, activation energy and reaction order. The 
initial estimates of these paramet^s were taken from the linear 



Asymptotic rote of cokingxlO^ kg/(m^)(h) 



Fig. 11 Determination of rote constant for asymptotic rate 
of coking at different temperatures 


Ink 
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Table 5: Effect of Process Variables on the Rate of 
Coking of Thiophene Free Naphtha 


Run 

No. 

T, K 

r ,3 

c Kg steam 
^ Kg naphtha 

lo4 Kg mole 

A 3 

m 

r xlo^ 

1 

1073 

0.20 

0.47 

1.39 

0.18 

2 

1073 

0,30 

0.47 

2.48 

0.33 

3 

1073 

0.40 

0.47 

3.07 

0.44 

4 

1073 

0.49 

0.47 

4.04 

0.53 

5 

1083 

0.20 

0.47 

2.47 

0.40 

6 

1083 

0.30 

0.47 

2.97 

0.49 

7 

1083 

0.40 

0.47 

3.39 

0.59 

8 

1083 

0.49 

0.47 

4.95 

0.99 

9 

1093 

0.20 

0.47 

3,61 

0,90 

10 

1093 

0.30 

0.47 

5,01 

1.23 

11 

1093 

0.40 

0.47 

4.88 

1.32 

12 

1093 

0.49 

0.47 

5.53 

1.50 

13 

1103 

0.20 

0.47 

5.06 

1.60 

14 

1103 

0.30 

0.47 

4.88 

1.65 

15 

1103 

0.40 

0.47 

5.88 

2.17 

16 

1103 

0.49 

0.47 

6.75 

2.35 

17 

1083 

0.20 

0.70 

1.79 

0,28 

18 

1083 

0.40 

0.70 

2.03 

0.32 

19 

1083 

0.49 

0.70 

'2,22 

0.36 

20 

1083 

0.49 

0,47 

7.37 

1.55 

21 

1083 

0.49 

0.47 

9.01 

2,00 

22 

1083 

0.49 

0.47 

9.83 

2.25 

23 

1083 

0,49 

0.47 - 


2.60 
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analysis. The optimization method used was a combined 
Gauss-Newton and quasi-Newton subroutine. The calculations 
were done on a DEC 1090 system and the detailed computer listing 
is given in Appendix C. Due to the strong correlation between 
the frequency factor and activation energy# a reparameterization 
was necessary. For the reparameterization# the rate constant 


was expressed as 


^o = 


E,1 1__, 

- R^T * T ' 
e m 


where T^ is the average temperature for all the runs. The 
value of kjj E and n obtained by nonlinear optimization together 
with the 95% confidence limits are shown in Table 6, The 
comparison between the experimental and calculated rate using 
this model is shewn in Fig, 13, The average percent deviation, 
calculated as 

) X 100 

was 6,0% and the maximum error between the predicted and 
experimental coking rate was 22%, Usirg non-linear regression 
the activation energy and reaction order were found to be 
58,5 kcal/mole and 1,32, respectively. Compared to the values 
reported by Kumar and Kunzru (1985), the activation energy is 
higher whereas the value of the reaction order is lower. The 
difference could be due to two reasons. Firstly, the metal 
surface used in the two studies were different and secondly, 
in their study, the water used for steam generation contained 


i z 

N ^ ^ 


^^a^cal, " ^^a^expt. 


^^a^expt » 
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Table 6: Parameter Estimates for the 
Non-linear Coking Model 


Frequency factor 

Activation energy 


Kq coke 

kc a 1/mole 

n RSS 

(„2)(hr).(£a_E2iS)n 



(1.53 + 0.96)xl0^^ 

58.5 + 7.5 

1.32 + 0.09 1.34x10*^ 
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50 ppm CS2, The presence of s in feed water can definitely 
influence the order of the reaction as was found in this 
study as discussed below, 

4*6 #2 Modell ing of Coke Formation During Pyrolysis of Thiophene 
Cont ainl nq Naphtha ; 

All the runs for rate of coking in the presence of thiophene 
were taken at a constant temperature. It was found that for 
thiophene— containing feed, the rate of coking could also be 
related to the aromatic concentration. The order of the 
reaction for coking was found by a linear plot of In r^ vs In 
as shown in Fig, 14, The order of reaction was found to be 
2,2, which show that the presence of thiophene increases the 
reaction order for coking. This order is comparable to the 
value of 1,98 reported by Kumar and Kunzru (1985) for CS 2 
present in the system. It seems that the radicals produced 
due to thiophene decomposition take part in the homogeneous 
and/or heterogeneous reactions to increase the order. Due to 
the large number of free radicals reactions taking place during 
naphtha pyrolysis, nothing can be said about the exact mechanism 
by which this occurs. 



Effect of inlet thiophene concentration on reaction 
order for rate of coke formation 




CHAPTER 5 


CONCLUSIONS AND RECQMMENDATT ON 


5,1 Conclusion : 

Based on the results of this study, the following 
conclusions can be made: 

1. The rate of coke formation is lowest on Inconel 600 
surface in comparison with SS 304 and SS 316. 

2, The presence of metal oxides formed on the 
reactor surface during decoking promote the 
rate of coke formation. 

3, Rate of coke formation Increases with an 
increase in aromatic concentration of the 
feedstock. 

4, Rate of coke formation increases with an increase 
In the reactor temperature and space time but 
decreases as the inlet steam to naphtha ratio 

Is Increased. 

5. Rates of coke formation decrease with increasing 
thiophene concentration of the feed. 

6. The coking model developed in this study in the 
absence of thiophene in the feed can be represented 
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’• -aptotic coU„, rate 

in the presence ot thiophene can he expressed 

approximately second order reaction 
Involving the aromatics. 

g-ecommendat-tnn - 

The foUc«l,^ recommendations are proposed for further 
Study* 


1. A detailed analysis of the and liquid 
products should be attempted to further 
improve the model, 

2. The naphtha pyrolysis can be conducted on 
swaged coils and the validity of the main 
pyrolysis and coking models can be checked 
for such coils. Moreover^ the effect of 
varying diameter may be studied for coking 
rates. 

3. Sulfur effect on naphtha should be checked 
more rigorously by analysing in detail the 
liquid and gaseous products and measuring rates 
of Coke formation at various temperature for 
naphtha mixed with thiophene. 
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APPENDIX A 


Calibration Details of Various Columns 


SI 

No. 

Component 

Peak area 
per C.C., 
(Duropak 
column) 

Peak area 
per c.c, 
(Poropak-Q 
column) 

Peak area per 
microlitre , 
(25% DECS on 
chr omo sorb -P 
CO lumn) 

1. 

Hydrogen 

- 

106666 

- 

2. 

Methane 


213333 

- 

3. 

Ethane 

- 

24166 

- 

4. 

Ethylene 

- 

15000 

- 

5. 

Propylene 

44,000 

- 

mm 

6. 

n**Butane 

79,500 

- 

- 

7. 

Benzene 


- 

107500 

8. 

Toluene 

■mm 

- 

117500 

9. 

Xylene 

- 

mt 

231610 


9 
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APPENDIX B 


Rates of Coke Formation with Run Time at Different 
T«nperatures and Space Time 


1, Temperature = 1073 K 


6 


n Ai Kg steam 
• Kg naphtha 


Space 
time, s 


0.20 


0,30 


0.40 


0,49 


Run time, 
h 


0.5 

1.0 

2.0 

2.5 


Rate of coking, X 10^ 

hr 

0.27 

0.23 

0,20 

0,18 


0.5 

1.0 

2.0 

2.5 


0,40 

0.36 

0,33 

0,32 


0.5 

1.0 

2.0 

2.5 


0.48 

0.48 

0,44 

0.44 


0,5 

1.0 

2.0 


0.69 

0.58 

0.50 

0,53 


2.5 
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2. Temperature *= 1083 K 


6 » 0.47 

Kg naphtha 


Space 
time, s 


Run time, 
h 


Rate of coking x 10'' 

la 


hr 


0.20 


0.30 


0.40 


0,50 

0,81 

1.0 

0.72 

2.0 

0.45 

2.5 

0,40 

0,5 

1.00 

1.0 

0,82 

2.0 

0,50 

2.5 

0.49 

0.5 

1.32 

1.0 

1,05 

2.0 

0.61 

2.5 

0.59 

0.5 

1.55 

1.0 

1.20 

2.0 

1,00 

2.5 

0,99 


0.49 
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3. T«nperature = 1093 K 

° • Kg naphtha 


2 

Run time. Rate of coking x 10 , 

h Kg 

hr 



0.5 

1.05 


1.0 

0,98 

0,20 

2.0 

0.92 


2.5 

0,90 


0.5 

1.50 


1.0 

1.37 

0,30 

2.0 

1.25 


2.5 

1.23 


0.5 

1.0 

0,40 2.0' 

2.5 

1.0 
2.0 
2.5 


1.70 

1.45 

1,34 

1.32 

2.10 

1.85 

1.52 

1.50 


Space 
time, s 


0.49 



60 



Space 
time, s 

Run time, 
h 

2 

Rate of coking x 10 , 

Kg 

m^ hr 


0.5 

1.80 


1.0 

1.75 

0.20 

2.0 

1.59 


2.5 

1.60 


0,5 

2.00 


1.0 

1.80 

0,30 

2,0 

1.68 


2.5 

1.65 


0,5 

2.65 


1.0 

2.40 

0.40 

2.0 

2.20 


2.5 

2,17 


0.5 

2.76 


1.0 

2.58 

0.49 

2.0 

2.33 


2.5 

• 

2.35 



APPENDIX C 


Computer Program for the Nonlinear Least 
Square Estimation of Model Parameters 



^ c » U=l 
)■ f i , 3 

i' ( n # r 

i ( I ''O = 1 


:b,u6 7f'>04 
= 1 ./ii03, 

, aS i:I"l >2 


tos, n=i •<.U,o,-*i 

ffis5lO^Ko2/ 

?‘ci^ii=b:7U?:-54 

Y(i 6 ) = 2 . ^ 51^702 

r(i7,i J=i./i08i, 

l'C17,2) = 1.7'3fc>y^ 

I. l?',n= •^^2 ^-04 

Y( i q 5 =ti , ibb-U? 
T(aDrU = l./108i. 
H20,2)s7.i7fc:-04 

iao5=1.55fc:-02 
rC 2 liu = l./l08^;.. 

rui# 2 ) = 9. 005^-04 

lcn5=2.ooK:S2, 

TC22l2)=9.B3E-04 

n225s2.2§e702 

T523j25 = i5,,07E-54 
SfC235=2.60E-02 


>‘8iTEi"2KS<?E^SU’'HEf' 


TEf iOuT# 999 nr ^ . .... 

RlTECNOUTjr99w63CXCJ)i^®lfp7;.^ 

ALL LSFUNHM,NfX,FVEe,FJAC,23) 

mol li. '*9 Ii 4 ■** '4H € ... M 



rE12.5) 

V 'UT^ri fft I* \ w *»• r - - - - ■'” * 

COfJTiMUE ^ « 

DO 121 Izht 

A?l!j)4ri'rJ)+F^^C(K,I)*FJACCK,J) 

:bK,t|.iiuE 


J=1 rN3 


0 ri 6 0 C) L ^ n r ^ ^ f 1 I "V • 1 5S 1 # f 

“.UtE(i'i(^UT,9994) 


u.4l' EHPOR EXII TyPE,l3,22rt ... 

^ "IK-Inn/Z/r^U! ON EXIT, THE SU'4 OF SQUARES IS/F12.4J 
t;l.y4HfCl|H &r the; FniMl',3Fl?.4) ^ ...i-i/.,/ 

K;n<'lAH6A4, 1 A3) 

F'lK'IATCI UX,&(2X,E12,5)) 
i '*' 1 J 


SEE RQiJTXMr. i.inCfh‘'.E'- , I H i’ ) 


)2 

5 


SihiROUriNK 
iMTEGER f.JC 


LSFON2CM,N,X,FVEC:,FJACCrLJC) 
M , N 


KFAF. FJACC(f.j5,NJ .fVLCCCM) ,X(N) 

REAL T(23,2),yC23) 
iNTFOKK i 
CO*1MOM YfT 
00 20 1 = 1 , M 
C=l,/ 1 D 0 B. 

FVeCCU)=X(U»EXP(XC2)*(r(l,l)-C))^(TCI,2)*^'X(i))-X( 1) 

FJACC(I,n = EXPCXC2)4'(TtI,l 3--C))*(T(I,2)^*XC3n ^ , 

FJACCa,2)=X(l)*CTCl,U-C)»EXP(XC23*(T(r,l}-C))*(T(I,2)**Xt31) 

FJACCCll3) = X(l )>^EXP(X(h*CT(I,l)-c5)»CT(i,2)*»X(3)l*AI,ar,(TCT,2J) 

CONTINUE 

DO 779 T = 1 ,M 

CONTINUE 

RETURN 

END 


TouoJSTng sub rou tin r*?! ^ s"*u u F’s r an d e rd "d e vIa F f q n 

^BROlffl NE~STAMV?x 7 FF 7 NrMrRSSQr 
DIMENSION X(3) ,FFCJ,33,FFTC3,3),FV'C3) 

INDIC=-1 
DO 10 1=1, N 
DO 10 J=1,N 
FFICI, J)=FFa,J) 

CAUL MATlN\/CNlFFT,Fy,INOIC,3) 

SI 3M Al=RSSD/( FLOAT (M-N)) 

DO 15 I=|,N 
DO 15 J=1,N 
IFCl.Eg,j5G0TD 20 
soio iS 

FrrCI,J)aSIGMAl=»'FFT(I,J) 

SD2aFFfCl.J) 

•Ri?i? 2253 olSl{xJlifs^i,RSSO 

CONTINUE 
RETURN 
END 


5.8,3X,'SDl=',Ei5.8,3X,'PSS=',El5.8J 


raBROUflNE"FOR IMVERSI0N’’0F MATRIX 


subroutine MATINVCM, AfX/INDlC,NRC) 

DIMENSION lROi<}(3 3 ,JCUI.C3) , JORD(3 3,XC3),ACNRC,MRC3 ,X(3) 
EPS=1,0E-15 

^'’AX=M 

if*'('ULE,20)r,OTU 5 
con 75 
D''l’KR=l. 

.)'} ! l< U-l,ri 

K..’l=K-t 



1 M 

ip-tr 
If-'C I. 


'■Vs* J J I I '•! 

TftC(U‘' = 1 ,KM1 


11 

11 


n 


JSCAMsj ,f(Ml 
E>^.lRQ«(TSCAN))GnTD 
, . Tp. JCOt.CJSCAN))GnTO 
w n 'I T I U E 

jsoTO ii 

PrvnrssACl.J) 

IR3^(K)=I 
j:3t,(K) = J 
CDMTINUE 

ICCABSCPI^DT) .GT.EPS)GDTO 
GO I'D 75 
IRa^KsIPGWCK) 

JCOLKsiCOnCK) 

DErF’HsPErKR^^PTVnT 
DO 14 J=1,MAX 

A(lROWK,J)=A(lROWK,j) /PIVOT 
CONTINUE 

A(£ROWK,JCUr,K) = I ./PIVOT 
DO IB 1=1, N 
AtJCK=ACI,JCOLK) 

IFCr.EQ.IROWKlGOTn 18 
ACi-JCOf.K)=-ATJCK/PlVnX 
DO 17 J=1,MRX 

IFCJ.ME, JCQLK)A{I,J)=AIT,J) 

CONTINUE 
CONTINUE 
DO 20 1=1, N 
IROWIalROWCl) 

JC0LI=JCQLC1) 

JORntIRO»^IJ=JCOLl , , 

IF(INDIC.GE,0)X(JCaLI)=R(IRO*iI,MAX) 

continue 

[NTCHaO 

im=N-i 

>0 22 1=1, NMl 

lP1=I+1 

>0 22 J=IPl,N 

.F(JORDCJ3,GE,JORD(I))GOrO 22 
rrSMPsJORDCj) , 
f0RDCj5=J0R0(l3 
fnRDCI)=JTEMP 
NTCH=INrCH+l 
ONTIWUE 

rClNTCH/2*2.NE.INTCH) DETERs-DETER 
FdNniC.GE.OlGOTO 26 
:0T0 75 
to 28 J= 1 ,N 
27 1 = 1, N 
.RJ'd = IRO''ni) 
rCOLl = JCOG(l3 
' r JC3LI)=AtTRnWI, J) 

:0''/nNUE 


AIJCK*AC1R0WK, J) 



i i U T) 

J” ji,j = jr'jr.l ,T) 

Vf ,1 ''/J ''J)=ACT,jCnLJ) 

Cl 'iTI r;ilK 

')'1 J= 1 ,H 

‘UI,J)=Y(.U 

' 1 '1 ^5 i '** I 

wfi rr.(^7,^ii?)(A(l,J).j=i 
FnHMATt/iDX, 'MATHTX ITiVRRSE 
f'TRMAt’l. 1 'JX,r,(4X,t:ji).8)/) 
m'fivw 
SfiD 


IS 


CALLED'/) 



